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Abstract 

The structure of an aqueous solution of sodium chloride at a planar elec- 
trode is investigated by integral equation techniques. With the central force 
water model the aqueous electrolyte is modelled as a mixture of sodium and 
chloride ions and partially charged hydrogen and oxygen atoms interacting 
via effective spherically symmetric pair potentials. The correlation functions 
obtained from the Ornstein-Zernike equation with Reference-Hypernetted- 
Chain closure gives a good description of the bulk structure (e.g. hydrogen 
bonded water network, solvation shell). With the bulk information and the 
Wertheim-Lovett-Mou-Buff equation we have calculated the density profiles 
at the uncharged and charged electrode. The rather rigid ice-like water struc- 
ture found previously at the neutral surface strongly repels the ions. Steric 
interactions between the differently sized ions and the ice-like water structure 
dominates the ionic distribution near the electrode. This model electrolyte 
also responds differently to opposite charges on the electrode. We found the 
asymmetry in the differential capacitance curve entirely determined by the 
response of the interfacial water structure. 
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I. INTRODUCTION 



In a previous article we have described an integral equation method for calculating accu- 
rate bulk correlation functions for the central force (CF) water model and reported results 
about the water structure at an electrode. The central force water model 0, which rep- 
resents water as a stoichiometric mixture of partially charged hydrogen and oxygen atoms 
interacting via spherically symmetric pair potentials, allows a completely molecular descrip- 
tion of water within the framework of integral equation theory. The effective pair potentials 
include the information about the proper triangular arrangement of atoms within the water 
molecule, the molecular dipole moment and the formation of hydrogen bonds. This water 
model has been successfully employed and optimized in computer simulations The 
water structure obtained near the electrode has been interpreted in terms of an ice-like 
arrangement of water molecules in good agreement with several computer simulation 
and integral equation results In this paper we extend our previous investigations to an 
aqueous model electrolyte and consider the consequences of the interfacial water structure 
on the structure of the electrical double layer [0. 

Early theoretical investigations of the electrical double layer neglected the molecular na- 
ture of the solvent. The solvent was approximated as a dielectric continuum and solvent 
effects appear only in mediating the Coulombic interactions. A mean-field treatement of this 
primitive model electrolyte leads to the Gouy-Chapman theory for the double layer, which 
holds in the dilute regime where the diffuse part of the electrical double layer is dominated by 
the long ranged Coulomb interactions. The remaining inner layer adjacent to the electrode 
was assumed to be formed by solvent molecules and excluding nonspecifically adsorbed ions. 
Several phenomenological models for the water monolayer employing adjustable parameters 
for the specific water interaction with the electrode were proposed to account for the exper- 
imental data 0,^. Of particular interest is the pronounced asymmetry of the differential 
capacitance with respect to the potential of zero charge, which characterizes the interface 
between an aqueous electrolyte and a mercury electrode . 



A step further towards a microscopic description of the liquid side of the interface was 
achieved by employing integral equation methods based on the statistical mechanics of liq- 
uids [|ll| or on density functional arguments |[T^JT3[] . First of all, in the 'primitive model' of 
an electrolyte the ions were treated as hard spheres with a point charge in the center. Then 
the continuum solvent was replaced by molecular models. Even for the simplest model elec- 
trolyte of equal sized hard sphere ions and dipolar hard spheres a detailed fluid structure 
emerges throughout the interface as a consequence of both solvent and ionic interactions 
1^-|T6| . With respect to the tetrahedral charge distribution within a water molecule, Torrie 



et al. supplemented the point dipole model by a quadrupole tensor of symmetry and strength 
appropriate to reasonably represent the water molecule. As a result of the highly orienta- 
tional interactions between these water-like molecules an ice-like solvent structure is formed 



near a smooth, curved and planar surface IjT^ in a region, which is usually attributed to 
the 'inner layer' of the double layer structure. Adding an octupole moment to the model 



yields a net polarization of the interfacial solvent structure at the neutral surface |T8|. In 
spite of some difficulty in making the electrode completely planar, the interfacial structure 
obtained with these multipolar water models resembles computer simulation results for more 
elaborate water models 0. 
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These more sophisticated steric water models (e.g. ST2, TIP4P and SPC) have been 
employed mainly in computer simulations because of their complexity. Several authors have 
focused their attention on how the detailed atomic structure of the electrode affects the 
interfacial structure of pure water for instance through chemical bonding, surface corrugation 
I9| , |20[] and image interactions modelling a metal electrode The model parameters are 



based on quantum mechanical ab initio calculations or on thermodynamic data. Present day 
computer simulations are restricted, especially for an aqueous electrolyte near an electrode, 
to systems with only a few ions and a few hundred of water molecules representing a highly 
concentrated electrolyte. Serious problems arise also from limitations in computing time, 
because long simulation runs are necessary to observe the ions entering and leaving the 
interfacial region |2^ , |23[| , and special molecular dynamics techniques have been applied p4 



Therefore it is difficult to obtain statistically reliable results from computer simulations near 
an interface. However, computer simulations for the complex electrode-electrolyte interface 
are just at the beginning. Recently, first results of an X-ray diffraction experiment indicated 



some interesting water structure near the electrode as well [25 



In comparision with the simulations, integral equation methods are more versatile when 
considering less concentrated electrolytes or when applying an external electric field on the 
interfacial structure and varying the field strength. This is in particular important for 
calculating the differential capacitance. In this paper we use the Wertheim-Lovett-Mou- 
Buff (WLMB) integral equation for determining the density profiles at the uncharged and 
charged electrode. We take the central force model as an appropriate water model for our 
calculations. The most dramatic effect of the specific water structure established near a 
planar surface on the electrical double layer is to be expected for ions with small cationic 
and large anionic size. Such a candidate is sodium chloride dissolved in water, which is 
particularly important in many electrochemical and biological applications. 

In Sec. II. the model for the aqueous solution of sodium chloride and the method for 
calculating bulk correlation functions is described. Results for the structure of the ionic 
hydration shell and for ion pairing are presented and discussed. In the first part of Sec. 
III. the WLMB-equation for the density profiles of a mixture of charged particles in contact 
with a planar electrode is given. The equation requires particle-particle correlations, which 
are taken from the bulk fiuid. In the region near the electrode the water solvent shows the 
ice-like layering. We discuss the effect of the ice-like water structure on the distribution 
of the differently sized sodium and cloride ions. In addition, we study the influence of an 
external electric fleld on the interfacial structure and the implications on the differential 
capacitance. Conclusions are summarized in Sec. IV. 



II. THE BULK ELECTROLYTE 

A. Model and computational method. The model electrolyte used within this paper is 
designed to represent sodium chloride dissolved in water. As an appropriate water model for 
the integral equations we employ the central force (CP) water model in the revised version 
10]. The central force model describes water as a stoichiometric mixture of partially charged 
hydrogen and oxygen atoms [qh = 0.32983 e and qo = — 0.65966 e) interacting via pairwise 
additive central force potentials, which are designed to reproduce the formation of water 



3 



molecules and their tetrahedral coordination via hydrogen bonds. The detailed form of the 
pair potentials are given in reference [|l|. Implicitly, the water molecules formed by the 
central force pair potentials are flexible. In spite of its simplicity this water model is known 
to yield good results for the bulk structure and for the water structure at a hard wall 

The interaction potential between the ion and each atom of the water molecule consists 
of a Lennard- Jones pair potential added to the Coulomb interaction 



+ 4e 



Ca/3 



12 



(1) 



We have fitted, according to the procedure proposed by Pettitt et al. the Lennard 



Jones parameters to the Hartree-Fock binding energy of the ion-water dimer in the equi- 
librium configuration, which has been determined for some alkali- and halide-water dimers 
by Kistenmacher et al. |^ with quantum mechanical ab initio calculations. The Lennard- 
Jones parameters for the interaction between a sodium or chloride ion and the atoms of a 
CF-water molecule are listed in table |. This parametrization of the ion-water interaction 
potential is by no means unique. Another set of pair potentials for describing the interac- 
tion between a sodium or chloride ion and the CF water molecule were derived by Bopp 
et al. [Q. Their pair potentials are able to yield a quite satisfactory representation of the 
dimer energy surface obtained from the ab initio calculations over a wide range of ion- water 
seperation. But they imply some important inconsistencies in the asymptotic behaviour of 
the correlation functions and as a result some thermodynamic properties, for instance the 
pressure calculated according to the virial equation or the isothermal compressibility, are 
not defined. 

The pair potentials between like ions are of the Huggins-Mayer form 



Ur 



r 



exp[-r/po 



(2) 



with potential parameters taken from |^. The interaction between the sodium and the 
chloride ion was modelled in addition to the Coulomb interaction by a Lennard- Jones pair 
potential with parameters eNaci = 2.816 x 10~^'^erg and <7j\faci = 3.6 A as proposed by 
Smith and Haymet. |2^. Taking instead the Huggins-Mayer form of reference for the 
pair interaction between the unlike ions and considering an ion concentration of 0.01 molar, 
convergent solutions of the Ornstein-Zernike equation are not obtained with the iteration 
procedure used for calculating the bulk correlation functions. We have found that conver- 
gence is not achieved because of a thermodynamic instability of the electrolyte solution, 
which can be characterized as a demixing of ions and the water solvent. The demixing 
cannot be avoided by changing the form of the interaction potential or the approximations 
employed in the iteration scheme, only the onset can be shifted [ pO[| . 

Since within the chosen model electrolyte the aqueous solution of sodium chloride is 
simply a homogeneous mixture of charged particles interacting via additional short ranged, 
spherically symmetric pair interactions, the Ornstein-Zernike equation for calculating the 
correlation functions reads 



Kp{r) 



(3) 
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together with the closure relation 



Ing^pir) = --^Uapir) + Kp{r) - Cap{r) - B^p{r). (4) 

kB is the Boltzmann constant, T is the temperature and denotes the number density of 
particle species 7. Including the so called bridgefunction Bapir) in the closure equation (|[) 
gives an exact relation between the pair potential and the pair distribution function 
Qap = 1 + hap- These bridgefunctions are known only for simple liquids and mixtures but 
not for more complex fluids like water or the aqueous electrolyte. Therefore these functions 
have to be approximated. Bap = is the familiar hypernetted chain (HNC) approximation. 
Application of the HNC-approximation for calculating the water correlation functions failed 
to reproduce the structural properties of pure water Previously we have developed 



bridgefunctions for the central force water model quite similar to those proposed by Ichiye et 
al. [^. Employing our bridgefunctions in the Ornstein-Zernike formalism yields a faithfuU 
description of the structural and thermodynamic properties of pure water under standard 
conditions. For instance two hydrogen atoms are bound to an oxygen atom forming together 
the triangular water molecule. Approximately two more hydrogen atoms are located near 
the oxygen atom at a distance, which is comparable to the hydrogen bond length indicating 
the formation of the hydrogen bonded water network. We have added these bridgefunctions 
in the closure relation (^) for the water correlations in the electrolyte, while the remain- 
ing ion-water and ion-ion correlation functions are calculated with the HNC-approximation. 
The set of coupled Ornstein-Zernike integral equations (^) together with the closure relations 
(I) are solved iteratively at fixed temperature T, water density and ion concentration. The 
results are the correlation functions hap{r) and Cap{r) in the bulk electrolyte. For details of 
the numerical calculation the reader is refered to reference In the calculations, presented 
in this article, the temperature is fixed at T = 300 K and the number densities for hydrogen 

o _3 o _3 

and oxygen atoms are chosen pH = 0.06690 A respectively po = 0.03345^4 equivalent 
to the mass density p = 1 gjcrn? . 



B. Structural results. With the integral equation method presented above we have calcu- 
lated bulk correlation functions for a 0.01 molar aqueous solution of sodium chloride, which 
means 5555 water molcules per sodium chloride molecule. In Fig.|T]a and |^ the pair distri- 
bution functions ga^ir^ of hydrogen and oxygen atoms around a central Na^- and C/~-ion, 
respectively, are shown. Both ion-oxygen distribution functions exhibit a sharp first max- 
imum determining the position of water molecules in the first hydration shell around the 
ion. We found the first oxygen atoms at a distance rj<iaO = 2.34A apart from the sodium 
ion and at rcio = 3.52A for the chloride ion. The positions of hydrogen and oxygen atoms 
in the ionic hydration shells are summarized in table |I|. 

Due to the stronger ion-water interaction the hydration shell of the smaller sodium ion is 
more pronounced. The broad maximum of the function gNauif) is obviously a superposition 
of two seperate maxima with different weight. We have fitted two gaussian functions to this 
part of the curve to determine the position of the two maxima. The first maximum is found 
at TNaH = 2.69A and the second at r^aH = 3.2lA. The hydrogen atoms contributing to 
the first hydration shell of the A^a"'"-ion are located at distances slightly smaller and larger 
than expected for a single water dipole perfectly aligned in the radial electric field of the 
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sodium ion. The arrangement of atoms described by the pair distribution functions can be 
interpreted as the water dipole in the iVa "'"-hydration shell tilted with respect to the oxygen- 
ion vector and one of the 'lone pair'-valencies of the water molecule is orientated towards the 
ion. Obviously, this orientation is energetically more favourable because one more hydrogen 
bond with water molecules in the second hydration shell can be established compared to 
the configuration where the water dipole is radialy pointing away from the sodium ion and 
always both 'lone pair'-valencies are directed towards the ion. Due to the packing of water 
molecules in the first hydration shell some molecules are also tilted with respect to the plane 
formed by the dipole vector and the oxygen-ion position vector. This leads to the double 
maximum of QNani^)- In order to evaluate the mean orientation of water molecules in the 
cationic hydration shell and to compare with computer simulation results it is convenient to 
consider the angle 6 enclosed by the dipole vector of the water molecule and the oxygen-ion 
position vector. We deduce the mean value of 6 from the positions of the maxima of the pair 
distribution functions gapir) under the assumption that the geometry of the water molecule 
in the hydration shell is not strongly perturbed by the neighbouring ion. This assumption 
seems to be well justified by computer simulations [Q. A value of ^ = 100° corresponds to 
the first maximum of gNaHif) and 9 = 145° to the second. With the weight of the gaussian 
fit functions we obtain 6 = 117° for the mean orientation of water molecules in the cationic 
hydration shell. Our results for 6 are summarized in table |T| together with results obtained 



in computer simulations. With regard to the results of Heinzinger et al. [28|, 9 is obtained 



from their values of {cos 6). The second set of values for the mean orientation is deduced 
from the pair distribution functions reported by Smith et al. [|^. Some uncertainties should 
be attributed to the results obtained from the pair distribution functions, since the angle 
6 is very sensitive to the geometrical parameters of the water model and to the location of 
pair distribution maxima. Nethertheless, 6 calculated with the integral equation method is 
in fairly good agreement with the computer simulation results. 

In addition we have calculated the running coordination numbers 

rr 

Naf3{r) = Att / dxx^gafsix) (5) 
Jo 

of hydrogen and oxygen atoms around the central ion in order to characterize the hydration 
shell. In Fig.|l]b and 0b the coordination numbers are plotted as a function of distance from 
the center of the ion. As a consequence of the pronounced A^a"''-hydration shell the function 
NNaoii^) possesses a well defined plateau. The coordination numbers, hsted in table ||, are 
defined as the value of Napir) at the first minimum of ^'^/^(r). Around the A^a+-ion we found 
on average 6.5 oxygen atoms. In our calculation the ratio N^^aH /N^ao of the coordination 
numbers according to the first minima of the pair distribution functions is 2.28, quite similar 
to a simulation result |33], indicating that also hydrogen atoms of the second hydration shell 
contribute to the first (7Ara//('")-iiiaximum. 

In a dilute solution the energy of ionic solvation consists of two parts: the solvent re- 
organization energy due to the creation of a polarized cavity in the pure solvent and the 
interaction energy between the ion and the surrounding water molecules already properly 
reordered. In the present work the ion-solvent contribution to the solvation energy is calcu- 
lated by determining the potential energy of the ion in its hydration shell according to 

POO 

[/-'- = 47r PP drr'gap{r)uap{r). (6) 

I3=0,H 
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The results for the Na'^- and the C/~-ion are presented in table is dominated 

by the Coulomb potential close to the ion. The orientation of water dipoles around the 
smaller A^a^-ion is more pronounced due to the strong electric field and therefore 11^°'''" 
is more negative for the A^a^-ion. Our value = — 8.60eV^ is close to the result of 

Chandrasekhar et al. ||3^ for the rigid TIPS-water model. 

The total work required to bring an additional ion into the system is given by the chemical 
potential fia- The excess part due to particle interactions can be determined from the 



two-particle correlation functions [|36| and in case of the HNC-approximation this relation 
reads 

f^a = 2vr J2 PP (^^^"^ [^lisir] - 9apir) Capir) - c«/3(r)J . (7) 

l3=H,0,ion ° 

In addition to the potential energy gain f/^°'^ of the ion due to ion-water interactions the 
excess part fi"^ of the chemical potential includes also the work required for creating a cavity 
for the ion in the water network and therefore must be less negative than 11^°''^. Our 
value = —3.23eV follows this prediction and is in good agreement with the value 

— 3.93eV^ determined experimentally for the free energy of solvation of a iVa"'"-ion 



Around the C/~-ion (cf. Fig.^a) hydrogen atoms are arranged in two distinctly seper- 
ated shells and the corresponding oxygen atoms are located in between. This distribution 
of atoms can be interpreted in terms of a water configuration in the anionic hydration shell 
where one hydrogen atom of the water molecule is directed into the liquid while the second 
is pointing towards the C/~-ion. Similar structural results were observed in analogous com- 



puter simulation studies |28,29| and in neutron diffraction experiments |38|. Due to such 



almost linear orientation of one oxygen-hydrogen bond towards the chloride ion each water 
molecule in the first C/~-hydration shell should contribute only with one hydrogen and one 
oxygen atom to the first gciuij)- and (yfc;o(?^)-maximum. This is expressed by the calcu- 
lated coordination numbers Ncih = 18.2 and Ncio = 17.9, although in comparision with 
the computer simulation result they are too large. This discrepancy may be attributed 
to the chosen ion-water pair interactions and to the HNC-approximation employed in the 
integral equation procedure. The coordination numbers are very sensitive on the details of 
the pair distribution functions, especially on the position and the magnitude of both the 
first (7Q,/3(r)-maximum and -minimum. As the size of the ion increases the coordination num- 
bers increases too because of an increasing volume available for water molecules in the first 
hydration shell. On the other hand the resolution of the first hydration shell decreases, as 
exemplified by the weak coordination plateau in C/~-hydration shell presented in Fig.^, be- 
cause of the weaker electrostatic attraction. As a result of the overpopulation of the anionic 
hydration shell, we obtain rather large values of f/^°- = — 7.34eV^ and /i^^- = —6.25eV 
compared with computer simulation results (cf. table [TVD and the experimental value of 
— 3.53eV^ for the free energy of solvation ||3^. With the RISM integral equation method 
Pettitt et al. |2^ calculated for the C/~-ion a more negative ion-water contribution to the 
solvation energy than for the A^a"'"-ion. This is clearly in disagreement with several com- 
puter simulations and contradicts the idea, that the solvation shell of the smaller cation 
should be bound more strongly because of the larger electrostatic attraction. The too short 
chloride-hydrogen distance might be the reason for their very large value of U^i-. 

The calculated structural properties of the ionic hydration shells for a dilute aqueous 
solution of sodium chloride are summarized in table |I| and compared with available experi- 
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mental data obtained by x-ray or neutron diffraction techniques and with computer 
simulation [28,^| and other integral equation results [26,31 ,^|. The results for the position 



of the first ion-oxygen and ion-hydrogen peaks are in good agreement with the experimental 
data and with those obtained by Smith et al. using the same CF-water model but different 
ion-water pair potentials. In the calculation of Thuraisingham et al. all bulk correlation 
functions, in particular for CF-water, are obtained with the HNC-approximation. This ap- 
proximation does not reproduce the structural features of bulk water. As a result the number 
of oxygen atoms in the first hydration shell is overestimated. An improvement is achieved 
in the present work by adding appropriate bridgefunctions in the closure relation for the 
water correlation functions. Also some results of an integral equation theory (RISM) for 
interaction site fluid models are added in table |T| ||2^ . 

Since in a dilute ionic solution only a few water molecules are interacting directly with 
the ion, the pair distribution functions QhH) Ooh and goo are indistinguishable from those 
of pure water reported in |T]]. 

The averaged structural effects of the solvent on the ion-ion interactions may be discussed 
in terms of the potential of mean force (PMF), defined as 

W^fsir) = -ksT \ng^fs{r). (8) 

Our results for the potential of mean force (cf. Fig.^ and Fig.^ show some characteristic 
structure due to the molecular nature of the solvent. In comparision with the bare pair 
interaction potential Ua/3{r), plotted in Fig.|^, the first minimum of WNaciif^) is drastically 
reduced in magnitude and slightly shifted towards larger r-values. This strongly reduced 
attraction between the Na~^- and CZ~-ion at contact arises primarily from a partial overlap 
of the hydration shells of the single ions. Debye-Hiickel screening is of minor importance. 
The first minimum of Wi\iaci{f) at r = 3. 42 A represents the 'contact ion pair' (CIP), whereas 
the second minimum at r = 5.36A corresponds to the 'solvent seperated ion pair' (SSIP), 
suggesting a solvation bridge between the sodium and the chloride ion. A combination of the 
configurations depicted in Fig.[l|a and Fig.|^ for one bridging water molecule explains the 
position of the SSIP-minimum. For the intervening water molecule a position is energetically 
more favourable, where the orientation matches those in the hydration shell of a single Na~^- 
and C/~-ion. Therefore the water molecule is slightly displaced from the line connecting both 
ion centers and the hydrogen atoms are orientated towards the CZ~-ion. This interpretation 



is in line with a previous molecular dynamics simulation , where the solvent configuration 
due to the minima of WNaCiii^) is inspected. 

The relative depth of the minima determines the stability of the CIP- and the SSIP-state. 
In our calculation with the CF-water model the first minimum of W]sfaci{T^) is found to be 
deeper than the second. The barrier for ionic dissociation (CIP —>■ SSIP) is 0.97 kcal/mol 
and for ionic association (SSIP — CIP) 0.30 kcal/mol. The barrier for the dissociation 
process is higher than for the association process and the CIP state is the most stable 
configuration. An indication whether the associated or the dissociated state is prefered by 
the sodium-chloride ion pair is provided by the ionic association constant [H^ 

WNaClir " 



Ka = 4:71 Na / dr exp 
Jo 



ksT 



(9) 



where denotes the Avogadro constant and R is the position of the barrier maximum 
seperating the SSIP- from the CIP-state. For the 0.01 molar solution of sodium chloride we 
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obtain Ka = 10.13 //mo/. In a dilute ionic solution with concentration c the coordination 
number Ka x c is interpreted as the fraction 6 of unlike CIP-pairs in the system. Recent 
MD simulations suggest the tendency for ion pairing to increase when the dipole moment 
of the water molecule decreases ^l|. Compared to the dipole moments of the water models 



employed in these simulations the mean dipole moment of the flexible CF-water molecule 
in the fluid phase ((/i) = 2.154 which is deduced from our pair distribution functions 
calculated with the HNC+B-approximation, is smaller. Our larger value for Ka, when 
compared with the simulation results, is therefore consistent with the higher probability of 
the CIP configuration when water models with a smaller dipole moment are employed. 

Despite the purely repulsive nature of the pair interaction between like ions, the respec- 
tive potentials of mean force, as plotted in Fig.|, show some oscillatory structure. The exis- 
tence of minima in the PMFs indicates ion pairing resulting from a subtle balance between 
the interionic repulsion and an effective attraction mediated by the ionic solvation shells. 
The PMF for Na'^ — Na~^ and for Cl~ — Cl~ are significantly different. The tendency for 
solvent-bridged sodium ion pairing is rather weak as indicated by the shallow minimum of 
WNaNa{f) at r = 4.52A and around 6A. These positions are in fairly good agreement with 



a previous computer simulation using the CF-model |2^. The first minimum of WNaNa{r) 
corresponds to a configuration where a water molecule is in a bridge position between the 
two sodium ions with the hydrogen atoms pointing into the liquid ESl. In contrast to the 



W]yaNa{r) OUT Calculated Wcicii^) possesses a prominent stable minimum at r = 6.6A. With 
the characteristic bond lengths in the chloride-water dimer we found this position consistent 
with a configuration where a water molecule is placed between the chloride ions orientating 
each of its hydrogen atoms towards one of the ions. The water molecule serves as a bridge 
in the chloride ion pair. The relative stability of the solvent-bridged chloride ion pair is 
known from neutron scattering experiments on highly concentrated electrolyte solutions. 
The Cl^ — Cl^ pair distribution function extracted from the measured partial structure 
factors shows a prominent peak at r = 6.lA for the aqueous solution of Ni2Cl and at 



6AA for the LiC/-solution [^. In addition to the 'dilute peak' around 6 A, a second 



maximum of gcici{f) is observed in neutron scattering experiments at r = 3.75A only for 
high ionic concentration. This peak is reminiscent of the molten salt regime and dimishes as 
the solution becomes more dilute. The RISM integral equation approach predicts only one 
stable minimum of the Cl~ — Cl~ PMF at r = 3.5A even at infinite dilution and obviously 
fails to reproduce the 'dilute peak' p6| . Although the positions of stable minima of the PMF 



are solely determined by the specific ion-water interactions, we found that the depth of the 
minima and therefore the stability of the ion pairs depend strongly on the ion concentration 
and the approximations employed in the integral equation calculations. 



III. THE ELECTROLYTE AT THE ELECTRODE 

A. Computational method. Several integral equations for calculating the density profiles 
at an electrode-electrolyte interface exist. A first relation between the particle density Pa.{^) 
and the potential energy V^(l) of a particle at coordinate 1 due to an external potential. 



e.g. the electrode, is given by the Born-Green- Y von (BGY) equation [EH] 
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Vilnp„(l) = - 



ViK(l) + ^ / d2p^(2)(7„^(l,2;[p])ViM„^(l,2) 



(10) 



where Uaisil, 2) is the pair potential and gaisi^, 2; [p]) denotes the pair distribution function 
in the inhomogeneous density distribution in front of the electrode. The functional depen- 
dence of ga/3 on the whole set of one-particle densities {p^(l)} is indicated by the notation 
da/si- [p])- By density functional arguments |1T2|JT3[| an alternative to the BGY-equation can 
be derived, where the density profile is related to the direct correlation function Cap instead 
of the total correlation function hap 



In 



Pa 



^K.(l) + E / ^2M2)-p^] £dXc^pil,2;[p,]). (11) 



This equation is exact too, but requires a functional integration of the direct correlation 
function over a sequence of density distributions [px] = {p^(l;A)} from the initial state 
(A = 0), e.g. the homogeneous fluid, to the final state (A = 1), when the external potential 
is completely switched on. By the same density functional arguments an even simpler 
relation between the density gradient and the external potential can be obtained, which 
avoids the functional integration. Originally this equation was derived by Wertheim and 



independently by Lovett, Mou and Buff |^8| for simple fluids, hereafter called the WLMB- 
equation. In the generalisation for a particle mixture the WLMB-equation reads 



Vi lnp„(l) = - ^ ViK.(l) + E / ^2 c„^(l, 2; [p]) V2Pp{2). 



(12) 



In contrast to eq.(|Tl|) the BGY- and WLMB-equation require only the inhomogeneous 
particle-particle correlations in the final density distribution near the electrode and not 
for a sequence [p\] of densities. Calculating the inhomogeneous particle-particle correlation 
functions from the Ornstein-Zernike equation requires great computational effort. Only for 
much simpler systems like the primitive model electrolyte in contact with a hard wall, this 
complicated calculation has been done For this reason several approximations for the 
inhomogeneous particle-particle correlations have been applied in the density profile equa- 
tions ([T^)-(|T^. Our experience has lead us to the conclusion that an approximation of the 
direct correlation function Cq,/3(1, 2; [p]) is more robust than of the total correlation function 
/ia^(l,2; [p]) p]| , P| . In the simplest approximation the inhomogeneous direct correlation 



function is replaced by its uniform-fiuid counterpart 

Ca/3(1, 2; [p]) ^ c„^(l, 2; {p^}) = c„^(l, 2). 



(13) 



Going a step further, the inhomogeneous direct correlation functions can be locally approx- 
imated by correlation functions of the homogeneous fiuid evaluated at appropriate weighted 
densities [p!2| , p0| -|5^ . For our complex model electrolyte we have achieved very good results 
using the approximation (|1^). Especially the ice-like water layers in the vincinity of the 
electrode were formed [0]. On this stage of approximation eq.(0) and the WLMB-equation 
(|I^) become equivalent. For numerical application we prefer the WLMB-equation. 

Considering a planar electrode perpendicular to the z-axis of a Cartesian coordinate 
system, the one-particle density Pa(l) and the external potential energy Vq,(1) are functions 
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of z only. The gradient is a derivative with respect to z and the integration perpendicular 
to the z-axis can be carried out seperately. This leads to the WLMB-equation in the form 

dlnpa{zi) 1 dVaizi) /-oo dpf3{z2 



2n r -^T^ r (14) 

o J— oo CLZo J\z-I—Z7.\ 



dz\ kj^T dz\ p J— oo dz2 J\zi—Z2 

Since our model electrolyte consists of charged particles, eq.(]T^) contains diverging terms 
due to the Coulombic interactions. For an analytical treatment of these long-range terms 
seperating the short-range part and the Coloumbic tail of the direct correlation function, as 
proposed by Ng [^, is appropriate 

/ N SRf ^ lo^ll^ erf{\r) 
Cap{r) = c^p (r) - , (15) 

where erf{Xr) is the error-function. With the definition of the charge density q{z) = 
J2aQoiPa{z) the diverging Coulombic terms are collected under the sum over /?. Requir- 
ing bulk neutrality {pa{z) ^=°° pa) and Pa{z < 0) = in the left halfspace, the charge 
density vanishes on the boundaries of the integral and integration by parts yields 



dlnpaizi) 1 dVa{zi) dpp{z2 



dZi ruB-t dZi ^ J— oo dZ2 J\z\—Z2\ 

27T_q, 



+ / dz2q{z2)sgn{zi- Z2)erf{X\zi- Z2\). (16) 

KrI J— oo 



In order to investigate the influence of an external electric field on the structural properties 
of the electrode-electrolyte interface, the elctrode is covered with a homogeneous charge 
density u: 



V^{zr) = V^''{zi)-27Tq^ujz,. (17) 

The repulsive nature of the electrode or some specific chemical bonding is incorporated in 
the short ranged interaction potential V^^{zi). Under the constraint of exact compensation 
of the surface charge density by charges in the fluid, the mean electric field due to the charge 
distribution in the right halfspace {zi > 0) 

POO 

E{zi) = -47r / dz2q{z2). (18) 

Jzi 

satisfies the boundary condition E{zi = 0) = Attuj at the electrode. With this assumption, 
using eq.(|T^) and eq.(|TB]) and introducing the complementary error-function erfc{x) = 1 — 
erf{x), equation ([T6|) can be rewritten in the form 



dlnpaizi) 1 dVa{zi) qaE{zi) foo dpfs{z2 



dZi KbI dZi KbI p J-oo dZ2 J\zi~Z2\ 



2rcqa 
kBT 



/oo 
dz2 q{z2) sgn{zi - Z2) erfc{X\zi - Z2\). (19) 
-00 



The information about particle-particle correlations are contained in the short ranged part 
of the direct correlation functions, calculated seperately from the Ornstein-Zernike equation 
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for the homogeneous mixture as discussed in section II. Starting with bulk density pa at 
^1 — s> oo an integration yields the density profile 



The set of coupled WLMB-equations ( ]T9| ) was solved iteratively applying a standard mix- 
ing scheme to stabilize the iteration procedure. In each iteration step the condition of 
charge neutrality has been carefully controlled. With this method the density profiles for 
an aqueous solution of sodium chloride in contact with an uncharged and charged elec- 
trode are calculated. The electrode was assumed unpolarizable. The short-range interaction 
V^^{z) = A/z^ {A = 1.88 X 10~^'^erg) of the electrode is purely repulsive and equal for all 
particles. This fast decaying repulsive potential is very similar to a hard wall, but conver- 
gence of the iteration scheme is achieved more rapidly for the 'soft' wall. 

B. Results for the uncharged electrode. Fig.^ shows the density profiles for sodium 
and chloride ions as well as for the hydrogen and oxygen atoms calculated for the 0.01 
molar solution at the uncharged, repulsive electrode. In front of the electrode there is a 
densely packed water layer. About 3A apart a second less densely packed layer is found. 
This distribution of hydrogen and oxygen atoms resembles the density profiles obtained for 
pure water at the same type of wall, which have been interpreted in terms of an ice-like 
water structure 0]. In comparision to these previous calculations the density of the first 
water layer is slightly reduced. Both ion species are strongly repelled by this densely packed 
water layer. In the region of the second water layer sodium ions are found with almost 
the same probability as in the bulk, whereas the larger chloride ions are still hindered from 
penetrating into this part of the interfacial region. Since the electrode acts equally repulsive 
on both ion species, steric and electrostatic interactions between the differently sized ions 
and the ice-like water matrix are responseable for the difference in the ionic density profiles. 

Information about the water structure is obtained by inspecting the charge density q{z) = 
q[pH{z)—2po{z)]. In Fig.^D this difference of particle densities is plotted for the water solvent 
in the electrolyte and for pure water. Both charge density profiles are very similar; that 
means the water structure at the electrode is not much perturbed by the presence of ions. 
The sequence of positive and negative charge density is characteristic for the ice structure. 
An schematic illustration of the ice structure can be found for example in reference [Q. Each 
water molecule provides two hydrogen atoms and two 'lone pair '-valencies in tetrahedral 
geometry. Since reorientation of water molecules without disrupting the hydrogen bonded 
network is possible under the constraint that each hydrogen bond is occupied by only one 
hydrogen atom, the ice structure has some internal polarizability. Thus the ice-like water 
structure can persist also under the infiuence of external electric fields. The charge density 
profile of water at a planar electrode has been already disscused in detail previously 
Here, we briefiy review the most important structural features: i) the broad density peak 
at contact actually consists of two oxygen layers with hydrogen atoms in between. In this 
double layer water molecules are connected by hydrogen bonds strongly tilted with respect 
to the surface normal. This yields a very small distance between the oxygen layers, ii) 
The detailed structure of another water double layer around 4A, where the negative charges 
dominates, is not clearly resolved in the charge density profiles, iii) The low density region 




(20) 
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between the double layers is characterized by the pronounced double maximum of positive 
charge density. Hydrogen atoms are arranged in two layers, forming hydrogen bonds between 
successive water double layers, iv) The calculated water structure at the electrode starts 
with a layer of oxygen atoms. Since our model electrode acts purely repulsive on all particles, 
it is energetically more favourable for water molecules immediately at the electrode to place 
their hydrogen atoms into hydrogen bonds within the first water double layer. Therefore 
the orientation of the first water layer is fixed by the water structure and this yields a 
spontaneous net polarization or surface potential of the electrode-water interface, even at 
the uncharged electrode. Due to the strong inclination of hydrogen bonds within the water 
double layer, the dipole moment of water molecules at contact is only slightly orientated 
towards the bulk fluid. 

The structural and electrostatic properties of the ice-like water significantly affects the 
distribution of the solute in the electrode-electrolyte interface (cf. Fig.^). It is of particular 
importance for the interfacial ionic distribution how the differently sized ions fit into the 
ice-like water structure. For instance, the smaller A^a"'"-ion (crjva+ = 2.4 A) penetrates into 
the second water double layer around 4 A. The hexagonal water coordination within a 
double layer is wide and flexible enough, so that the ion is able to slip into the region 
between the first and second water double layer. In this region the water structure has 
still enough flexibility to match the geometry of the iVa"'"-solvation shell without disrupting 
the hydrogen bonded network. The seperation between the first and second water double 
layer increases slightly. Simultaneously positively charged hydrogens are transfered from 
this interlayer region into both neighbouring water double layers (cf Fig.pD). But the ice- 
like sequence of positive and negative charge density still persists. The redistribution of 
positive charges can be interpreted as a reorientation of water molecules to optimize their 
interaction with the neighbouring sodium ion. Water molecules of the upper and lower 
double layer locally provide an environment for the A^a"'"-ion, which resembles the bulk 
solvation shell and is also compatible with the ice structure. Therefore the local density 
PNa{z) is close to its bulk value. In the region between the first and the second water 
double layer the mean electrostatic potential (cf. Fig.^ becomes repulsive for the A^a"'"-ion 
due to the large number of positively charged hydrogen atoms arranged in hydrogen bonds 
connecting successive water double layers. Although the iVa+-ion, when further approaching 
the electrode, is attracted by the net dipol moment of the first water double layer, penetration 
into the densely packed first water double layer can be only accomplished by a strong, local 
distortion of the ice-like ordering. But in this region the ice structure is rather rigid and the 
approach of the ion is therefore strongly hindered. 

The behaviour of the larger C/~-ion {aci~ = 3.24 A) in the electrode-electrolyte interface 
is quite different. The chloride ion is too large for the hexagonal rings of hydrogen bonded 
water molecules, already in the partially disordered second water double layer. In order 
to provide a proper solvation environment for the Cl~ -ion hydrogen bonds have to be dis- 
rupted, which is energetically unfavourable and therefore the density pci{z) decreases in the 
interfacial region. Although the mean electrostatic potential (cf. Fig.^ becomes attractive 
for the C/~-ions in the interlayer region, the rigidity of the ice-like water structure and thus 
the steric repulsion increases. These competitive interactions yields the local maximum of 
Pci{z) around 4 A. Additional chemical ion-surface or image interactions with a metal elec- 
trode can overcome the repulsion of the rigid water layer [^. At a more realistic electrode 
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contact adsorption is expected for the larger CI -ion 
The mean electrostatic potential across the interface 



/■oo roo 
^(Zl)=/ dZ2E{z2) = ATI dZ2{Zi~ Z2)q{z2) (21) 



is plotted in Fig.^ for the electrolyte and for pure water. The shape and the magnitude of 
both curves are very similar throughout the entire interface. Adding ions to the water reduces 
the potential drop across the first water double layer, indicating a decreased surface dipole 
moment due to the reorientation of water molecules in order to optimize their interactions 
with sodium ions in the neighbourhood. The contribution of the ionic charges to '^{z) is 
rather small because of the very low ion densities, especially near the electrode. The total 
dipole moment due to the charge density in the interfacial region is given by the overall 
electrostatic potential drop across the interface 

i.e. the work required to move a test charge from the bulk liquid with '^{z = +00) = 
to the electrode. In the literature the potential drop \l/o at zero surface charge density 
is sometimes denoted as the 'potential of zero charge'. In electrochemistry the measured 
'potential of zero charge' (PZC) means a different quantity (cf. Fig.|^). Often the 'surface 
potential' x is used, which is related to \l/o for the uncharged electrode by the sign conven- 
tion X = ~^o- Our value \I'o = —TAmV for the electrode-electrolyte interface is somewhat 
smaller in magnitude than \E'o = — 110m\^ calculated previously for pure water [Q. The 
negative sign of \l/o reflects the preference of the water oxygens close to the surface. Our 
surface potential for water at a planar, inert wall agrees in sign with the surface potential 
for the liquid- vapor interface predicted from theoretical and experimental data [Q, and 
with the surface potential obtained by Torrie and coworkers for a water-like solvent model 
with multipolar hard spheres [ITS. Recently Booth et al. have calculated \l/o according to 



eq. (p2[) using eq. ([TI|) with bulk correlation functions for determining the density profiles of 
hydrogen and oxygen atoms for the central force water in contact with a quasi-hard wall 
||55| . Their value \l/o = -l- 600 mV has the opposite sign. From the corresponding curve 
for the electrostatic potential, presented in Fig. 6a of reference together with the Pois- 



son equation, which relates the curvature of the mean electrostatic potential to the charge 
density q{z), it follows a reversion of the surface dipole moment in disagreement with our 



and other results The reason for this reverted orientation is not quite clear, because 
Booth and coworkers used a model and computational treatment completely equivalent to 
our procedure. Only the particle-particle correlations taken from the bulk might be differ- 
ent. According to our understanding, the driving force for the spontaneous polarization of 
the electrode-water interface is the tendency of the interfacial water structure to maximize 
the number of hydrogen bonds |l|]. 



C. Results for the charged electrode. In this section we discuss the influence of an ex- 
ternal electric field due to an electrode covered with a constant charge density on the distri- 
bution of ions in the electrode-electrolyte interface. The surface charge density u is increased 
up to = ±10A fiC/cm^ corresponding to an electric field strength of 1.18 x 10^ V/cm, 
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comparable to a very strong electric field applicable in electrochemical experiments. Our 
previous calculations for pure water have shown that the ice-like water structure near the 
electrode is preserved under the influence of such field strengths. The reorientation and 
the exchange of hydrogen atoms between the two positions in a hydrogen bond provides 
for enough polarizability. In case of a small negative surface charge density alignment of 
water dipoles with respect to the electric field is hindered. For the hydrogen atoms of the 
contact water layer it is energetically more favourable to remain still between the oxygen 
layers as long as the energy gain in the external potential is smaller than the binding energy 
in the hydrogen bond. Therefore the polarizability decreases. Only when the hydrogen 
bond energy is overcompensated by the potential energy gain of the hydrogen atom in the 
external field, some hydrogens follow the attractive electric field. As a result the interfacial 
polarizability increases. This reorientation of water molecules starts at u = — 3.7 nC/cm^ 
Experimental and theoretical investigations found a maximum of the interfacial entropy 
for surface charge densities in the range from —4 to —QpC/cm? |]^, where the hydrogen 
atoms of the contact water layer, according to our calculations, can occupy two different 
positions. The mechanism of water reorientation still holds for the water solvent in the 
electrode-electrolyte interface. The ion density profiles, presented in Figure have to be 
discussed with respect to this specific behaviour of water molecules in an applied electric 
field. For a surface charge density > the electric field cannot force the large C/~-ion 
to penetrate into the contact water layer, whereas an electric field of opposite sign but of 
same magnitude brings the A^a+-ion close to the electrode surface. This different ability of 
approaching the electrode is a result of two competing forces: the attraction due to the ex- 
ternal electric field and the steric repulsion of the differently sized ions by the ice-like water 
matrix established in front of the electrode. When the external field is attractive enough for 
the hydrogen atoms of contact water molecules to flip towards the electrode, the A^a"'"-ions 
finally invade the first water double layer and even get into contact with the surface (cf. 
dashed-dotted line in Fig.^). A measure for the excess of ions in the electrode-electrolyte 
interface is 



r 



a 



/ dz[pa{z) - Pa]. (23) 

Jo 



The plot of Fq, for the A^a"''-ions vs —u in the inset of Fig.^ja reflects the rapidly increasing 
penetration of iVa"'"-ions into the surface water structure beyond u ~ —ApC/cm'^, where 
the pronounced water reorientation becomes significant. Contact adsorption of A^a"'"-ions 
has been also observed in computer simulations [^. 

In order to elucidate the effect of the asymmetric response of ionic charges in an applied 
electric field on the differential capacitance C = duo / d'^ q{uj) , we evaluated '^o{oj) across the 
interface for a series of surface charge densities. Some of these values \l/o are summarized 
in table 0. For comparision the last column of table shows the values of \I'o(^^) for 
pure water. As already mentioned, at zero surface charge the electrostatic potential drop 
across the interface is only slightly smaller than for pure water, but the negative sign is 
retained. The same trend was observed by Torrie and coworkers when adding ions to their 



multipolar, water-like solvent model []T8|. In the third column we have listed the change 
A\l/o = ^'o('^) — ^I^o(^ = 0). Only for small surface charge densities |A\l/o| behaves like 
a symmetric function of cj, whereas for large negative surfaces charges |A^o| grows more 
strongly than for corresponding > 0. This implies a smaller differential capacitance for 
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the negatively charged electrode. In order to decide whether this asymmetry is due to 
the response of the water solvent or due to the ions, we have examined the water and 
ionic contributions to the mean electrostatic potential \l/o(i^) seperately, as shown in table 
^ For the 0.01 molar solution the ionic contribution \l/jomc is rather small in comparision 
with "^water- FoT almost all surface charge densities "^ionic is independent of u. Only the 
pronounced adsorption of A^a"'"-ions yields a distinct increase in "^ionic- Since "^ionic is rather 
small and flat, the asymmetric response of the total "^oi^) is primarily a property of the 
interfacial water structure. Nethertheless, to some extent water is indirectly affected by the 
presence of ions through a restructuring of the contact double water layer, as emphasized 
by comparing the solvent component water and \E'o(OM) for pure water. water for the 
electrol5d;e varies less strongly with respect to u. As a result the differential capacitance 
for the electrode-electrolyte interface is larger than for pure water over the whole range of 
considered surface charge densities. 

The differential capacitance C = duj / d'^ q{uj) obtained from the derivative of a cubic 
spline fit to the data is shown in Fig.|^ as a function of surface charge density uj. The 

curves for the 0.01 molar aqueous solution of sodium chloride and for pure water, taken from 

are very similar. Because d'^ ionici!^) / doj is negligible in comparision with d"^ water{^^) / duj , 
the total capacitance of the electrode-electrolyte interface is dominated by the response of 
the interfacial water structure. The minimum in the differential capacitance curve has been 
identified with the beginning of water reorientation in the adsorbed water layer [0]. The 
asymmetry in the differential capacitance curve with respect to zero surface charge and the 
existence of a minimum at negative surface charge density are typical of the experimental 
data for simple electrolytes |T0|. It has been argued that the asymmetry is a consequence of 



ionic polarizability or of ion-specific interaction with the electrode [56|. Our result suggest 



that the minimum can be explained by the asymmetric water response near an electrode, 
at least for low concentrated electrolytes. We found the differential capacitance for the 
electrolyte larger than for pure water over the whole range of surface charges considered, 
indicating an overall larger polarizability of the interfacial structure. From the rather small 
direct effect of the ion distribution on the differential capacitance, we infer that due to the 
presence of sodium and cloride ions in the interface water molecules are already reorientated 
and less strongly fixed within the ice-like water structure yielding a larger polarizability of 
the interfacial water structure. 



IV. CONCLUSIONS 

We have extended our previously developed integral equation method for calculating bulk 
correlation functions and density profiles for the central force water model to an detailed 
microscopic model of an aqueous electrolyte. Spherical symmetric pair potentials for the 
additional ion-ion and ion-water interactions are designed to represent the different size of 
sodium and chloride ions as well as the orientation of water molecules around both ions. 
Only for the water correlation functions we have added to the HNC-approximation appro- 
priate bridgefunctions, which are optimized in order to yield the bulk water structure 
The calculated ion-water pair distribution functions reflects specific features of the ionic 
hydration shells in aqueous solutions. On average six water molecules are located in the 
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first A^a "'"-hydration shell and their dipole moment is tilted with respect to the ion-oxygen 
radial vector. The water molecules around the larger C/~-ion are bound more weakly and 
therefore the first hydration shell of the anion is less sharply defined. We have found that 
in the anionic hydration shell the water molecules are orientated with one of the hydrogen 
atoms towards the CZ~-ion. 

Taking these correlation functions of the homogeneous fluid for describing the particle- 
particle correlations near the electrode, we have calculated the density profiles for the various 
particle species near a planar electrode, which is neutral or charged and has some hard wall 
repulsion, equal for all particles. The charge density profile of the water solvent that is 
established at a neutral electrode resembles that for pure water, which has been interpreted 
in terms of an ice-like arrangement of water molecules [||]. This layered structure extends 
about two water double layers into the liquid. Water molecules adjacent to the surface are 
slightly orientated toward the bulk fluid and thus the interfacial water structure is polarized, 
even without any external electric field. This specific water structure induces an asymmetric 
ion distribution in front of the electrode depending strongly on the ion size. The smaller 
A^a"'"-ion can easily penetrate into the second water double layer and even comes into contact 
with the surface for sufficiently attractive fields. Penetration of the larger C/~-ion into this 
interfacial region is strongly hindered. Strong local, steric interactions between the ions 
and the ice-like water structure dictates this different behaviour. Electrostatic effects are 
only of minor importance. The cavities formed by the puckered hexagonal water rings 
within a water double layer are wide and flexible enough to match the size of the smaller 
Na~^-ion. Water molecules are reorientated without disrupting the ice structure in order to 
provide an solvation environment for the A^a"'"-ion. This reorientation of water molecules is 
responseable for a reduction of the surface potential from \l/o = — 110 mV in case of pure 
water to \l/o = —74mV calculated for the 0.01 molar solution of sodium chloride at the 
uncharged electrode. The CZ~-ion is too large for the rigid hexagonal water cavities and 
therefore the C/~-ions are strongly repelled. 

Under the influence of an external electric field of reasonable strength the ice-like water 
structure near the electrode is still preserved. Chloride ions are less effectively attracted by 
a positively charged electrode than sodium ions by an equivalently attractive electric field 
because the rigid first water double layer resists the electric field. On the other hand Na~^- 
ions suddenly pentrate into the contact water layer, when the electrode charge is large enough 
{—AjjiC/crn?) to induce 'flipping' of adsorbed water molecules within the ice structure. 

Although sodium and chloride ions behaves differently at oppositely charged electrodes, 
this has a negligible effect on the interfacial differential capacitance C = duj / d'^ q{uj) . At 
least for dilute aqueous electrolytes the asymmetry in the differential capacitance curve 
with respect to zero surface charge density is solely determined by the ice-like water struc- 
ture as discussed more extensively in reference fl]]. Due to the presence of A^a"'"-ions some 
water molecules are less strongly fixed within the interfacial water structure and already 
reorientated in order to optimize their interaction with the neighbouring ion. As a result, 
we have obtained a larger polarizability and therefore a larger differential capacitance for 
the electrode-electrolyte interface compared with pure water. Our present model has to be 
further developed for instance by introducing specific electrode-ion interactions in order to 
understand more features of the experimental differential capacitance curves microscopically. 
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FIGURES 

FIG. 1. a) Pair distribution functions gNaHir) ( ) and QNaoi^) ( — ) for a 0.01 molar aqueous 

NaCl-solution. b) The corresponding coordination numbers NjyaHi'i') ( ) and N]\faoii') ( — )• 

FIG. 2. a) Pair distribution functions gcinir) ( ) and gcioi''") ( — ) for a 0.01 molar aqueous 

NaCl-solution. b) The corresponding coordination numbers Ncih{t) ( ) and Ncioii") ( — )• 

FIG. 3. a) The potential of mean force for the ion pair Na'^ — Cl^ and b) the bare interaction 
potential UNaCiir)- 

FIG. 4. The potential of mean force for the ion pairs Na'^ — Na'^ ( ) and Cl^ — Cl~ ( — ). 

FIG. 5. a) The density profile Pa{z)/ Pcx for hydrogen ( ) and oxygen atoms ( — ), for Na^- 

(• • •) and C/^-ions (- • -) at the uncharged electrode, b) The charge density q{z) = q[pH{z)—2po{z)] 
for the water solvent in the electrolyte ( ) and for pure water ( — ). 

FIG. 6. The mean electrostatic potential ^{z) near an uncharged electrode for pure water 
(-•-), the 0.01 molar aqueous solution of sodium chloride ( ), ionic ( ) and solvent contribu- 
tions (• • •). 

FIG. 7. Schematic diagram of contributions to the electrochemical potential of an electron. 

Different contributions are plotted as —e(j){z). Electrostatic potential ( ), potential in front of 

liquid surface ( ), ^q: the surface potential of the electrolyte eq.(p2|), PZC: potential of zero 

charge measured in electrochemistry, <I>j: workfunction of the electrodes, xm'- surface potential of 
metal, S: total surface potential at reference electrode, AP: absolute electrode potential at reference 
electrode related to redox reaction. 

FIG. 8. The density profiles for the Na'^- and CZ~-ions at the a) negatively and b) positively 

charged electrode for surface charge densities u; = ( ), ±1.6 ( — ), ±3.2 (• • •) and ±6.4/iC/cm^ 

(-•-). The inset shows Tjsia{^) as a function oi —uj. 

FIG. 9. The differential capacitance for the 0.01 molar aqueous NaCl-solution and for pure 
water. 
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TABLES 



TABLE II. Structural properties of the ionic hydration shell at T=300 K. 







Ch- 


This work 

Clmi + Vi M til R I9qI 

Pettitt et alj] ^ 
X-ray @ 
neutron diffraction 


ion-oxygen dist 
2.34 

9 '^'^ 
z.oo 

2.30 

2.38 - 2.40 


ance in A 

3.52 
3.45 
3.2 - 3.34 


Ttiis work 
Smith et al. Q 

X^eLLlLL cL ai. 

neutron diffraction 


ion-hydrogen d 
2.69 and 3.21 
3.0 
o.vo 


istance in A 

2.78 
2.53 
9 n 

2.2 - 2.26 


This work 
Bopp et al.| §8| 
Thuraisingham et alj^ Q 
fchiye et al.| jsij 
Pettitt et al. [||| 
X-ray Q 


coordination m 

6.5 

5.9 
15.9 

5.7 
4.28 

4.6 


imbcr ion-oxygen 

17.9 
8.4 
17.3 

11.73 


This work 
fchiye et al. ^ 
Pettitt et al. 
neutron diffraction 


coordination ni 
14.9 

6.97 


imber ion-hydrogen 

18.2 
14.8 

5.3 - 6.2 



"MD, infinite dilution 
''RISM, infinite dilution 
'^MD, 2.2 m NaCl 
"'HNC, 0.069M NaCl 
''RHNC, infinite dilution 
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TABLE I. Ion-water Lennard-Jones parameters, cr in A and e in 10 ^"^erg. 



a 




O'aH 


EaO = £aH 


ci- 


2.222 
3.378 


1.547 
2.703 


0.1853 

0.03708 



TABLE in. Orientation of water molecules, mean values of @ in degree. 





Na+ 


ci- 


This work 
Smith et al.g l2| 
Heinzinger et al.|^ |28| 
Heinzinger et al.|^ |28| 


117 
124 
139 
132 


44 
34 
55 

54 



"MD with CF- water, from ga/si''') 

^MD with CF-water, from = cos~^(cosB) 

'^same as b with the ST2-water 



TABLE IV. Energy of solvation in eV. 





TTSolv 


jjsolv 

^ci- 


This work 
ChandrasekhaiQ et al. 34] 
Zhu^ et al. m 
Pettittf] et al. [|26| 


-8.60 
-8.47 ± 0.04 
-6.92 ± 0.53 

-7.21 


-7.34 
-6.21 ± 0.04 
-4.51 ± 0.53 

-8.25 



''MD, 0.444 molal 
''MD, 1.791 molal 
'^RISM, infinite dilution 
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TABLE V. Surface potential ^'o, A^o = ^o(w) — *o('^ = 0), ionic and solvent contributions 
in volts for th(^ 0.01 M a(HUM)us NMC'l-soluliou and \['o(Oil/) for i^urc water. 



uj \iiCjcrr?\ 


^0 




^ water 


^ ionic 


^-olOM) 


-10.4 


-1.039 


-0.9652 


-1.018 


-0.0210 


-1.0890 


-6.4 


-0.670 


-0.5960 


-0.657 


-0.0130 


-0.7140 


-3.2 


-0.371 


-0.2973 


-0.359 


-0.0117 


-0.4119 


-2.4 


-0.296 


-0.2227 


-0.285 


-0.0116 


-0.3364 


-1.6 


-0.222 


-0.1483 


-0.210 


-0.0115 


-0.2608 


0.0 


-0.074 


0.0000 


-0.062 


-0.0114 


-0.1100 


1.6 


0.073 


0.1470 


0.085 


-0.0114 


0.0402 


2.4 


0.146 


0.2199 


0.158 


-0.0114 


0.1149 


3.2 


0.219 


0.2924 


0.230 


-0.0114 


0.1895 


6.4 


0.504 


0.5775 


0.516 


-0.0114 


0.4847 


10.4 


0.849 


0.9226 


0.860 


-0.0112 


0.8464 
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Electrode 



Electrolyte 



J I I I L 



10 



Reference 
Electrode 



